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LIST OF SYMBOLS 
A 2 

A cross-section area in m 

E q output voltage of a strain gauge bridge 

4 

I section moment of inertia in m 

4 

J section polar moment in inertia in in 

K gauge factor of a strain gauge, i e , unit 
resistance change per unit strain 

M bending moment m m-lb 

Q torque m m-lb 

R resistance m ohms 

T thrust load m lb 

V bridge excitation voltage 

e strain m microstrams, i e , micro-inches 
per inch 

Cg bending load induced strain 
£q torque load induced strain 
e.p thrust load induced strain 
a stress m psi 


u Poisson's Ratio 
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1 Introduction 

This project, initiated m September 1971, has as its objective 
the determination of propeller design parameters which yield minimum 
pro’peller noise, given constraints on propeller performance 

The study addresses both the theoretical and the experimental 
aspects of the problem Theoretical efforts have been directed at 
applying variational techniques to the noise equations Experimental 
efforts have been concerned with measurement of the propeller performance 1 
and noise of a Lockheed Y0-3A aircraft 

The following sections presents a report of the current state of 
this study 
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2 Theoretical Studies 

The objective of the theoretical aspect of the program is to deter- 
mine the optimum blade loading to minimize noise level, given the thrust, 
diameter, rotational speed and forward velocity. Both an aerodynamic model 
and an acoustic model for the propeller are required for this analysis 

2. 1 Aerodynamic Model 

The aerodynamic model at present consists of a vortex system util- 
izing a lifting line and a helical vortex sheet for each blade The blade- 
wise circulation distribution r (y) is the sought for quantity m the noise 
minimization problem 

A computer program has been written for the aerodynamic model based, 
essentially, on the analysis of Moriya (Ref. 21). The program, as it now stands 
assumes constant bladewise section aerodynamic characteristics. Modification 
to allow variable section aerodynamics is a relatively easy task and is now 
m progress The results of the aerodynamic analysis express the propeller 
performance m terms of the Fourier coefficients for the circulation distri- 
bution T(y) A draft of this analysis is enclosed as Appendix A. Some 
further shakedown and check out of the program is required before it becomes 
fully operational 

2. 2 Acoustic Model 

The acoustic model is a simplification of that represented by the 
Ffowcs -Williams -Hawkings equation (Ref 2.2, Eq 3-16) Presently only the 
force terms are retained, resulting m a model that is essentially Gutm's 
model with radially distributed line dipoles representing the blades 
Following successful treatment of this model, volume displacement effects 


will be included. 
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The sound intensity may be expressed m normalized form, as 


I* 


fi fi 

O O 


F ( r* , rp Y (t*) Y ( rp dr + dr^ 


where y is related to the circulation and F is a known function The con- 
straint of constant thrust is expressed as 
1 

y (r*) (l-r*)dr* = 1 
o 


Applying variational analysis to these equations provides a necessary condi- 
tion for I* to be a minimum 

f 1 h 

j F(r*,rpyCrpdrJ = - |r*(l-r*) 

o 

where b is a Lagrange multiplier. 

Initial attempts at numerical inversion "of this integral equation 
utilized both a Legendre Polynomial expansion and a Chebyschev polynomial 
expansion for y(r*), coupled with a Gaussian integration procedure, reducing 
the Integral equation to a matrix equation 

AL = C 

where L represents the unknown values of y at selected points on the blade 
Numerical inversion of this matrix equation was not successful In all cases 
tried, the coefficient matrix A was ill-conditioned 

The current approach involves a Fourier transform of the integral 
equation, implying the assumption that the noise is periodic m the blade 
frequency, with a Doppler shift This is consistent with the acoustic far- 
field assumption 

Some success is indicated with this approach, however convergence 
is exceedingly slow and, as a consequence, required computer time is large 
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Appendix B contains some details of the earlier analysis for a 
simplified case m idiich only the thrust component of blade lift was included. 

Appendix C contains some details of the current analysis utilizing 
Fourier Transforms. 

2 3 References 

2 1 Monya, T., "Selected Scientific and Technical Papers," University 
of Tokyo, Tokyo, 1959 

2.2 Goldstein, M E , "Aeroacoustics," National Aeronautics and Space 
Administration, Washington, D C , 1974 
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3. Experimental Activities 

All experimental efforts have been directed toward m-f light measure- 
ments of propeller torque and thrust and fly-by sound levels, utilizing a 
Y0-3A aircraft 

3 1 Aircraft Preparation 

Since the Y0-3A aircraft was received m disassembled, worn and 
damaged condition, a good deal of time and effort has been spent rendering 
it airworthy. Much of this effort was expended m obtaining maintenance 
and operating manuals, shop drawings, etc. The cooperation m these endeavors 
of Mr Harold Schuetz, USAAVSCOM, St Louis, and Mr. David Schnebly, Lockheed, 
Sunnyvale, is very gratefully acknowledged 

Major repairs have included a top overhaul on the engine, repair of 
a fractured vertical fm-fuselage attachment fitting, repair of aileron 
and wing tip damage, replacement of elevator attachment pins, manufacture 
of special jigs for checking drive shaft/prop shaft alignment and belt ten- 
sion, modification of the forward instrument panel, installation of VHF radio, 
removal and repair of the muffler system, and other items Photographs of 
the aircraft and the propeller drive systems are shown m Figures 3 1 and 
3 2 

3. 2 Measurement of Propeller Performance 
3 2.1 Introduction 

An essential element m the development of a low-noise propeller is 
the accurate measurement of its propulsive performance This means that the 
thrust available, the thrust power available, and the shaft power required 
must be measured, preferably m flight For the subject studies, a 
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Lockheed YO-3A aircraft is available for use as a test platform A signifi- 
cant part of the program has been devoted to instrumenting this aircraft 
for the m- flight measurement of thrust and power 

Two approaches to these measurements were considered m this study 

a) Airflow Measurement 

By pressure surveys of the flow field behind the propeller, 
determinations of the thrust and torque acting on the propeller 
are possible (Ref 3.1) However, this method requires correc- 
tion for the interference effects of the fuselage The accuracy 
of this method was judged insufficient for the purposes of this 
pro] ect 

b) Direct Thrust and Torque Measurements 

The Y0-3A's propeller is mounted on a cantilevered pro- 
peller shaft (Figures 3 2 and 3 3) This shaft is driven 
through a speed reduction drive system employing pulleys and 
belts Because of this drive system, the engine is isolated 
from the shaft and load cell measurements on the engine mounts 
would not yield propeller thrust data Such measurements could 
produce approximate torque measurements if pulley, belt and 
propeller shaft bearing losses could be determined 

As shown m Figure 3. 3, approximately 9-mches of the pro- 
peller shaft is easily accessible The strains induced m this 
shaft are direct measures of the thrust and torque transmitted 
to the propeller Thus, direct instrumentation of the propeller 
shaft (say, with strain gauges) can, m principle, yield the 
desired information This approach was selected 
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There is one major difficulty involved m this approach, 1 e , 
the thrust related strains are much smaller than the torque re- 
lated strains and are very small m an absolute sense No mechani- 
cal modifications to the shaft (e g weakening the shaft m ten- 
sion) were attempted because of airworthiness questions involved 
Instead, it was decided to attempt to measure the shaft strains 
directly 

322 Propeller Shaft Loads 

The Y0-3A aircraft is powered by a Continental IO-360A engine rated 
at 210 hp at 2800 rpm It drives a 3-bladed, constant -speed propeller 
through a sheave and belt system that provides a 3 33 1 speed reduction 
ratio. The blade pitch is controlled by a Woodward governor that supplies 
pressurized engine oil to the propeller through the hollow propeller shaft 
Based on data given m Reference 3 2 and assuming no drive system 
losses, it is estimated that a maximum of 186 hp is available at the pro- 
peller shaft to drive the propeller This, together with a propeller shaft 
speed of 841 rpm, yields a maximum torque to be transmitted by the shaft 
of 13,944 m -lbs. 

Again based on data from Reference 3 2, it appears that the thrust 
required m level flight, at sea-level, will be m the range from 250 to 
400 lbs. Although performance data are not available for the propeller, 
the assumption of an efficiency of 50% at the sea-level stalling speed of 
61 knots yields an estimated maximum thrust available of about 500 lbs. 

The thrust load instrumentation has been designed on the basis of a maximum 
thrust load of 500 lbs, but can measure thrust loads up to 1000 lbs 

Besides the torque and thrust loads, the shaft is subject to loads 
from three other sources 
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a} Bending Due to Propeller Weight 

Bending stresses due to the weight of the propeller are estimated 
to be due to a 1500 m -lb moment acting about the location of the 
instrumentation on the shaft 

b) Internal Oil Pressure 

The internal oil pressure produces an axial tension stress due 
to the action against the hub piston. It also produces a hoop 
stress (tension) that, through Poisson’s ratio, induces a com- 
pressive axial stress The net effect is a small axial tension stress. 
This net stress is about 28 psi for an oil pressure of 100 psi; this 
is equivalent to an additional thrust load of about 53 lbs 

c) Centrifugal Force 

At a shaft speed of 841 rpm, the outer surfaces of the 
shaft has an acceleration of about 24 g's The resulting axial 
stress is negligible compared to the other shaft stresses 
Table I summarizes the stresses and strains m the propeller shaft due to 
thrust, torque, and bending Both the centrifugal stresses and the net 
oil pressure stresses as calculated were negligible, however, their effects 
will be evaluated during calibration of the instrumentation The ratios 
of the maximum strains taken from Table I are 

£ T : e M £ 0 . 1 6 6 , 42 4 

max max rniax 

These ratios show the crux of the instrumentation problem, i e,, the very 
low thrust strains m the presence of much larger bending and torque strains 
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323 Strain Transducer Arrangements 

The sensor choice was dictated by the need to distinguish between the 
combined loads present, as well as the need to measure them reliably and 
accurately The sensor must also function over the range of environmental 
conditions met m aircraft flight. Strain gauges were chosen as the sensors 
Most strain gauge work is done m the 50-500 microstram range, 
within this range it is possible to measure changes of 1 or 2 microstrams 
However, the thrust strains on the shaft are 8. 85 microstrams full scale, 
see Table I. This low strain level poses problems m the thrust measure- 
ment The torque and bending strain levels are high enough to present no 
difficulties m reliable measurement by foil strain gauges 

Figure 3,4 outlines the decision process followed m choosing the 
method of measurement and selection of sensors Two mam avenues of approach 
were made for the thrust measurement 

1) Intensification of the low level thrust strains by mechanical 
means, using foil strain gauges as sensors. 

2) Direct shaft strain measurement of thrust strains Foil 
strain gauge and semi-conductor strain gauge evaluations were 
made and compared The semi-conductor gauges are 75 times as sen- 
sitive as foil gauges but have a high inherent temperature effect 
on sensitivity as compared to foil gauges 

Wheatstone bridges, with an active strain gauge m each arm were 
chosen If all four arms experience equal temperature changes, and if all 
four gauges are perfectly matched then temperature effects are automatically 
cancelled. Figures 3 5 and 3.6 illustrate the directions of the strains 
produced by the thrust, bending, and torque loads and 'the bridge arrangements 

used for measuring them. Each of the strains can be measured independently 
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of the other two and, m each case, the bridge output exceeds the output of 
a single active gauge 

When the shaft is subjected to thrust, bending and torque, the in- 
duced strains act along certain directions, as shown m Figure 3.5 Thrust 
and bending loads produce tensile and compressive strains m the axial 
direction Due to the Poisson effect, axial tensile and compressive 
strains generate compressive and tensile lateral strains. The shear strains 
induced by torque cause tensile/ compressive strains on mutually perpendicular 
axes inclined at 45° to the shaft axis, the sum of the components of torque 
strains m the axial and latenal directions are zero The sums of the 
thrust and bending strain components, m the torque strain directions, are 
not zero but their effects on the torque bridges will cancel 

Figure 3 6 shows the three arrangements of strain gauges used to 
measure thrust, bending and torque on the shaft, taking advantage of the 
directional properties just discussed Referring to this figure, we have 
Thrust. E q = [2(l+u)e T ] V 

Torque E q = [4e Q ] y 

Bending E q = [ 4 Eg ] v 

Ignoring second order effects, non-linearities and cross -sens\tivi- 
ties of the gauges, and assuming the gauges are perfectly matched and aligned 
m the required directions, each of these bridge arrangements reacts only 
to one type of loading and is insensitive to the other two loads 

3.2 4 Mechanical Intensification of Thrust Strain 

Larger thrust strains can be obtained by adding a load path parallel 
to the propellant shaft and incorporating m this added path a short, weak 


link 



Such an arrangement is shown m Figure 3. 7. The propeller shaft (1) carries 
most of the thrust load The parallel load path (made up of the support link 
[2) and the weak link (3)) carries only a small part of the thrust load, the 
amount is determined by the relative strengths of the weak link and the pro- 
peller shaft However, both load paths stretch equally under the load 
Because most of the stretching of the parallel path occurs in the short weak 
link, its strain is much greater than that of the shaft, l e , 

An analysis of this strain intensifier arrangement yields the following 


intensification ratio 

£, 


A,E £ 9 A E 

1 + V ( A^ + 0 


'-A P 
A 2 b 2 


where 


e = strain, micro-inches per inch 
£ = length 

A = cross-sectional area 
E = modulus of elasticity 

The subscripts 1, 2, and 3 refer to the components as shown m Figure 3.7 

and e^ is the shaft strain m the absence of the intensifier If the tensile 
o 

strength of the weak link is made much less than those of the support link and 
the propeller shaft such that 


and 


^1 ^ 3^3 

^ (rr) <K 1 
0 11 


£ 2 A 3^3 

(/) Cjrr) « i 

o A 2 b 2 


then the intensification ratio is approximately 


h 
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The system used to study the performance of the strain mtensifier is 
shown m Figure 3.8. The system consisted of two weak link assemblies, 180° 
apart, mounted on a full-scale model of the exposed portion of the propeller 
shaft. The links were fabricated of aluminum and were each instrumented 
with 4-arm strain gauge bridges. Using the smallest gauges available, a 
link length of 0.2 inches was achieved, yielding a length ratio, of 

40. The thickness and width of the links were 0.010 and 0.400 inches, 
respectively. The tensile strength ratios were 


A E 

X 2 — = 0.003 
l fc l 


A 3 E 3 

A 2 E 2 


= 0.001 


The test assembly was installed m a lathe (Figure 3.9) and could be 
loaded simultaneously m torque, bending, and thrust. Loadings up to 40% 
of full-scale were used. Typical results of the thrust calibration are 
shown m Figure 3.10, and the overall thrust and the torque interactions 
with thrust results are summarized in Table II. It is evident that the 
torque interaction is very large. Over the full-scale range of loads to be 
met m flight, the weak link system would have an output, due to torque 
interaction, of 860% and 212% of the full-scale thrust output. (The 
interaction results of the two links are different due to different 
arrangements of the gauges on each link. This was done to check the 
effect of gauge arrangement on the torque interaction.) Because of the 
large torque interaction effects and because the links required very close 
machining tolerances and were found to be very difficult to install and 
adjust satisfactorily, the mechanical strain mtensifier was rejected as 
a practical solution to the thrust measurement problem. 
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A second attempt at mechanical intensification was made using a deflec- 
tion sensor manufactured by DSC, Inc In this device, a small cantilevered 
arm is instrumented with semi-conductor strain gauges to measure the bending 
strains produced by deflection of the end of the arm The deflection sensi- 
tivity appeared very attractive for the present application, i.e , about 
6yV/V/y-mch compared to a shaft strain of 8.85 microstram for the 500 lb 
thrust load 

The sensor's arm has a flexure near its end to reduce the effects of 
bending moments other than those due to linear deflections It was hoped 
that this feature also would help reduce the torque and bending interactions 
m the thrust measurement Tests of several different installation arrange- 
ments showed, however, that the interactions were much worse with this sensor. 
Because of this, it was rejected as a possible thrust sensor 

3.2 5 Direct Shaft Strain Measurements 

With the abandonment of mechanical strain intensification, a study of 
the practical problems of direct thrust strain measurement was undertaken. 

For comparison purposes, both foil and semi-conductor strain gauges were used 
In the case of the foil gauges (which have low strain sensitivity but good 
temperature characteristics), both 120 ohm and 350 ohm gauges were used 
The latter gauges were used to take advantage of the higher excitation 
voltages possible. The semi-conductor gauges (which have high strain sensi- 
tivity but poor temperature characteristics) were 1000 ohm gauges with 
gauge factors of about 155 

The gauges were mounted on the model propeller shaft and were in- 
stalled m the lathe (see Figure 3 9) and loaded m combined thrust and 
torque. The torque sensor was made up of 350-ohm foil gauges fabricated 
as special purpose torque rosettes Typical results are presented m 
Figures 3.11 and 3.12 The results are summarized m Table II For 
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comparison purposes, the expected sensitivity of the foil bridge was 
E 

= .023 micro-volts/volt/lb thrust 

and the expected sensitivity of the semi-conductor bridge was 
E 

= 1 78 micro-vo lts/volt/ lb thrust 

From the results of Table IX it was concluded that the foil gauges 
are adequate for the torque measurement However, the higher sensitivity 
of the semi-conductor gauges was found necessary for reliable measurement 
of the small thrust strains Since the aircraft will produce a noisy en- 
vironment for the instrumentation system and since the sensor signals will 
be transmitted from the rotating shaft through slip rings, it is important 
that the sensor outputs be as large as possible to insure a high signal-to- 
noise ratio 

The selection of semi-conductor strain gauges for the thrust measure- 
ment was made m spite of their poor temperature sensitivity characteristics. 
The resistance of these gauges increases while the gauge factor decreases as 
the gauge temperature increases These temperature effects can be reduced 
by placing a suitable resistance (called a span resistor) m series with the 
bridge This resistor has a negligible resistance change with temperature 
As the bridge resistance increases (reducing the circuit current), the voltage 
drop across the span resistor decreases and voltage drop across the bridge - 
increases This increase m bridge (excitation) voltage tends to compensate 
for the decrease m bridge sensitivity The value of the span resistance 
can be used to control the temperature range over which compensation is 
achieved. 

The effectiveness of temperature compensation was studied by installing 
an electric resistance heater inside the model shaft and measuring bridge 
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output as a function of bridge temperature (which was monitored by thermistors 
mounted near the gauges) . Some typical results from this study are presented 
m Figure 3 13. These results show that the uncompensated bridge has a 
temperature sensitivity of 12.3 yV/V/°F That is (for a bridge initially 
balanced at 78 °F) at a bridge temperature ofl04°F, the temperature induced 
output would be about 40% of the full-scale thrust output However 
with a 915 ohm span resistor, this bridge, at 104°F, will have a temperature 
induced output of only 5% of the full-scale thrust output. In flight, the 
propeller shaft temperature will change due to changes in both ambient 
temperature and oil temperature inside the shaft. Thus, temperature com- 
pensation of the semi-conductor bridges will be mandatory In addition, sur- 
face temperatures will be measured so that corrections m bridge output can be 
made. 

The results presented m Figure 3 14 and Table II, show that there 
is a large torque interaction bridge output. This interaction possibly can 
be explained by one of the following arguments Theoretically there should 
be no torque strain components m the axial and lateral directions , which 
are the directions m which the thrust gauges are mounted This is true 
provided that the thrust gauges are exactly aligned on the shaft axes, and 
that the shaft is isotropic and homogenous In the present case, the ratio 
of the maximum thrust and torque strains is 1 42 4 Thus, small misalignments 
of the thrust gauges from the shaft's axial and lateral directions would cause 
a pickup of torque strain components that might be sizeable compared to the 
small thrust strains. The term misalignment here is meant to include both 
the deviation of the gauges from the axial and lateral directions caused m 
gauge mounting, as well as the deviation of the shaft axis from its theoretical 
direction. 

For the full-scale loads, the thrust and torque interactions per 



degree of gauge misalignment have been calculated to be* 


A maximum torque strain interaction of 96 8% 
of the maximum thrust strain would occur for 
each degree of misalignment of the thrust 
gauges 

A maximum thrust strain interaction of 2.23% 
of the maximum torque strain would occur for 
each degree of misalignment of the torque 
gauges 

Another possible source of this interaction could be the cross- 
sensitivity of the gauges This is normally a very small effect, but in the 
present case it may be significant because of the 1 42 4 ratio of maximum 
thrust and torque strains. 

For both the foil and semi-conductor thrust bridges, whose results 
are presented in Table II, the torque induced output was of a sign opposite 
to the thrust output. (The sense of the applied torque was the same as that 
which would occur m the aircraft ) When the sense of the torque was reversed, 
the sign of the torque induced strains also changed Both the foil and semi- 
conductor gauges were mounted on a shaft which had buckled slightly m an 
earlier test Before the gauges were installed, the shaft was remachmed, 
but the buckling could have caused a curved shaft axis 

Another test shaft was made and a new set of semi-conductor gauges 
were mounted on it When this shaft was tested under the same conditions as the 
older shaft,' the torque interaction m the thrust bridge output was still 
present but its maximum value had been reduced to 44 3%. Also, torques 
applied m the same sense as that m the aircraft, produced torque induced 
strains of the same sign as the thrust induced strains This result is 
opposite to that observed with the older shaft This result tends to support 
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misalignment as a cause of the interaction It is very difficult to exactly 
predict the misalignment effect because four individual gauges are involved 
However, for both the old and the new shafts, the torque interaction was 
calibratable and repeatable 

A bending interaction m the thrust output was also observed. This 
interaction was much smaller than the torque interaction The bending inter- 
action can also be caused by misalignment, as well as by small differences 
m the sensitivities of the gauges making up the thrust bridges For the 
new shaft, the maximum interaction was only 7 05%. This bending interaction 
also was found to be calibratable and repeatable 

3.2 6 Instrumentation Development and In-F light Systems Description 

3 2.6.1 Introduction and Background 

A block diagram of the thrust/torque measuring system is presented 
in Figure 3 15. The elements of this system can be grouped into three sub- 
systems, (I) strain measurement components mounted on the propeller shaft, 

(2) components located m the cockpit, and (3) auxiliary components. By 
way of background, a short discussion of the development of the strain 
measurement system is presented first 

Because of the very low thrust -induced strain levels, a high gam 
system with a high signal-to-noise ratio is required if the thrust measure- 
ment is to be successful Initially an AC system, with a 1000 Hz carrier 
frequency, was designed and built. This system exhibited sufficient sensi- 
tivity and was quite linear However, the zero-strain bridge output could 
not be nulled, even when the bridge was resistance balanced Furthermore, 
this zero load minimum bridge output level was sufficient to saturate the first 
amplification stage. Is is thought that the problem was due to phase un- 
balance between the semi-conductor bridge arms In view of this problem. 





3-14 


the system was redesigned to use DC excitation throughout The DC system, 
which will be used for m-flight testing, has been built and is now being 
checked out on a full scale model of the propeller shaft. With the exception 
of the data recorder, all components of the instrumentation system receive 
power from a regulated +_ 15 VDC power supply operating from the 400 Hz air- 
craft electrical system 

3.2 6.2 Components on the Propeller Shaft 

t 

The thrust, torque, and bending transducers (4 -arm strain gauge bridges) 
are mounted directly on the propeller shaft The associated electronics are 
mounted on a printed circuit board that is attached to the front face of the 
slip-ring assembly [see Figure 3 16 and Section 3263 below) Bridge exci- 
tation is by a regulated, IC 5 VDC power supply that is part of the shaft 
instrumentation package and operates from the 15 VDC source This additional 
regulation was selected to provide a stable source for the strain transducer 
circuits 

The bridge outputs are processed by IC differential amplifiers feeding 
into operational amplifiers [See Figure 3.17). The amplifier combinations 
are adjusted to give overall gams of 1000 The amplifiers are mounted m 
high thermal resistance packages and have active internal temperature con- 
trols giving very tight temperature control This arrangement provides 
highly stable gam and low temperature drift over a range of ambient tempera- 
tures from -55°C to + 125°C. 

3 2 6.3 Slip Ring Assembly 

A sectional view of the slip ring assembly is shown m Figure 3.18 and 
a photograph of the components, before final assembly, is presented m Figure 
3 19. The brush/ring material combination was selected on the basis of 
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information given m Reference 3 3 and after some preliminary laboratory 
checks The objective was a low resistance contact and very low slip ring 
noise The materials selected were hard brass rings and silver-graphite 
(50% silver by weight) brushes Two brushes and two internal leads are pro- 
vided for each ring to further improve the contact characteristics. 

In order to avoid the complications of a split slip ring assembly, 
the inside diameter of the assembly is 5 inches, large enough to slip over 
the flange at the forward end of the propeller shaft (see Figure 32) A 
split spacer is clamped to the shaft, to bring the diameter up to 5 inches 
, Two aircraft quality bearings separate the inner (rotating) and outer (fixed) 
assemblies. The outer assembly is thus carried by the shaft and is not 
rigidly attached to the aircraft frame. A soft rotational restraint will 
be used, minimizing vibration effects. 

3 2.6 4 Cockpit Components 

An amplifier control box and a data recorder are located m the 
observer's (forward) cockpit The control box contains active filters, 
having 3db roll-off at about 65 Hz, and output amplifiers for use m 
adjusting the overall gam of the system (See Figure 3 20) It will also 
contain auxiliary components discussed m Section 3 2.6 5 below 

The data recorder is a Gulton Industries, 8-channel, strip chart 
recorder provided by the University’s Aviation Research Laboratory (ARL) 
where it has been used for m-flight testing. It is powered by a 400 Hz 
to 60 Hz inverter, also supplied by ARL 

A bundle of shielded co-axial cables will connect the slip-ring 
brushes to the amplifier control box 

3.2 6 5 Auxiliary Instrumentation 


In order to correct the semi-conductor bridge output for temperature 
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effects, thermistors will be mounted on the shaft surface close to the strain 
gauges They will be connected to the cockpit control box through the 
slip rings and the cable Associated electronics will be contained on the 
printed circuit board on the shaft and m the control box Their output 
will also be recorded on the 8-channel recorder 

The apparent thrust load will contain a contribution due to the oil 
pressure m the shaft This pressure will be monitored by a pressure trans- 
ducer in the hydraulic line between the governor and the propeller shaft. 

Its associated electronics also will be m the control box with output 
recorded on the 8-channel recorder 

Finally the control box will contain a single stage timer to provide 
a time reference for the recorder output 

3.3 Fly-by Noise Measurements 

Capability of fly-by noise measurements has been established 
Measuring equipment consists of a General Radio 1962-9601 microphone with 
associated pre-amps and power supply, feeding into a NAGRA IV-SJ tape 
recorder. Analysis of the results is accomplished by use of a General Radio 
1564-A Sound and Vibration Analyz-er The noise recording and analysis 
equipment is University supplied 

Procedural checks of the noise measuring program have been satis- 
factorily completed. 

Weather conditions with sufficiently low ambient noise level to 
allow detection of the very low noise level of the Y0-3A aircraft have not 
prevailed at this writing This aspect of the program is ongoing and air- 
craft noise measurements will hopefully be achieved m the near future 
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Figure 3.1 The Y0-3A Aircraft with Cowl Removed 
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Figure 3.2 The Y0-3A's Propeller Shaft and Hub 
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Figure 3.4 Approaches to the In-F light Measurement of Propeller Performance 
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Figure 3.5 Propeller Shaft Strain Gauge Arrangements and Strain Directions 
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Figure 3.7 Schematic of Mechanical Strain Intensifier 




Intensifier I 



Intensifier II 



Figure 3.8 Mechanical Thrust Strain Intensifier, Mounted on Model Propeller Shaft 

3 


"Up" Component of Torque 
Load; Cable and Pulley 

d. 



"Down" Component 
of Torque; Cable 


Bending 
Load; Rod 


Thrust Load; 
Cable and 
Pulley 


Figure 3.9 The Combined Loading Arrangement for the Model Shaft 
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Figure 3 13 Temperature Compensation of the Semi-conductor Thrust Bridge 
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Figure 3 14 Semi-Conductor Thrust Rridpe Calibration for Tornue Tnfer'icfi on 



Figure 3 15 Block Diagram of In-Flight Instrumentation System 
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Figure 3.16 The Instrumentation System 
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Figure 3.17 Schematic of One Typical Channel of On-Shaft Instrumentation 
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Figure 3.18 Sectional Sketch of Slip Ring Assembly- 
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Figure 3.19 The Slip Ring Assembly 
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PROPELLER SHAFT LOADS 
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APPENDIX A 


NUMERICAL ANALYSIS OF PROPELLERS 


by R. M. Plencner 



LIST OF SYMBOLS 


axial interference factor 

rotational interference factor 

coefficients of sine senes 

aspect ratio 

number of blades 

non-dimensional chord, C/R 

chord, feet 

lift coefficient 

drag coefficient 

diameter, feet 

induction factor 

advance ration, V/nD 

loading factor 

rotational speed, rev/sec 

torque, ft-lbs 

torque coefficient 

distance along blade, feet 

radius of propeller, feet 

vector from origin to segment of trailing vortex filament 
vector from origin to a blade section 
thrust, lbs 
thrust coefficient 

velocity induced at the propeller disk, ft/sec 
forward velocity of propeller, ft/sec 



normal induced velocity- 

vector sum of rotational and forward velocities 
non-dimensional distance along blade, r/R 
rearward distance from the propeller, feet 

absolute geometric angle of attack 

induced angle of attack 

effective or sectional angle of attack 

blade angle 

circulation 

interval around the singularity 
efficiency 

parametric variable of the helix, radians 
tip speed ratio, V/J1R 
density, slugs/ft 3 
propeller solidity, BC/2irr 

angle between the relative velocity and the horizontal pi 
rotational speed, rad/sec 



I. INTRODUCTION 


A great deal of effort and expertise went into modeling the aero- 
dynamic characteristics of propellers m the early part of this century, 
However, with the coming of the jet age m the early- 1950's, the pro- 
peller was made obsolete for many applications and new theoretical work 
on propeller theories was almost non-existent Therefore, most of the 
work done m this area was undertaken before digital computers were de- 
veloped Since the problem of flow around a propeller is by nature very 
complex, these early approaches to the problem necessitated many simpli- 
fying assumptions and often made use of graphical techniques to obtain 
solutions to the problem 

Recent years have seen a renewed interest m propeller modeling for 
predicting their performance in a wide variety of applications These 
applications include more efficient propellers to reduce fuel consumption, 
propellers suitable for STOL aircraft, low noise propellers for commercial 
and general aviation aircraft, as well as rotors for wind mills. Many of 
the applications require more accurate prediction over a wider range of 
operating conditions than is possible with classical propeller theories 

The purpose of this paper is to review the classical approaches to 
the propeller problem and then dispense with as many simplifying assumptions 
as possible m order to produce a theory that more closely models the 
actual physical phenomenon 



II. CLASSICAL THEORIES 


Momentum Theory - A propeller produces thrust by increasing the 
velocity of a large quantity of air. This production of thrust is 
associated with a loss of energy which is due to the increase m kinetic 
energy, rotational motion imparted by the torque and frictional losses 
of the propeller blades In the simplest form of the momentum theory, 
first set forth by R E Froude fl] and W. J. Rankme [2], only the 
axial component of momentum is considered; the rotational and frictional 
losses are ignored. Therefore, the theory describes the characteristics 
of m ideal propeller. 

The axial momentum theory is developed [3, 4] by replacing the 
propeller with a propeller disk that has an infinite number of blades, 
and the thrust uniformly distributed over the disk. The flow is assumed 
to be incompressible and irrotational m front of and behind the propeller 
disk. Under these conditions the change m pressure across the propeller 
disk is equal to the total pressure head m front of the disk minus the 
total pressure head behind the disk. It can then be shown that at the 
propeller disk the increase m velocity over the free stream value is one 
half of the increase m the ultimate wake. The efficiency which is 
predicted by the momentum theory, which, is the maximum theoretical efficiency, 
is given by 


ri = propeller efficiency = 


CD 


where V is the free stream velocity and v is the velocity increase at the 
propeller disk. 



Practical applications of the momentum theory are limited due to 
the gross simplifications that are made The theory neglects rotational 
energy losses, non uniform thrust loading losses, profile drag, blade 
interference, and compressible drag changes Its direct application, 
therefore, is only principally useful m obtaining an upper limit on 
efficiency or a crude first estimate for the performance of a propeller. 
A more general form of the momentum theory [S] includes rotational 
momentum as well as axial momentum. However, this added complexity does 
not increase its usefulness m most applications. 

Blade Element Theory - The blade element theory [5] analyzes the 
aerodynamic characteristics of the propeller by estimating the aero- 
dynamic forces experienced by each blade element as it moves through 
the air. The blade element theory is evolved by assuming that the 
aerodynamic forces acting on a given blade element are equivalent to the 
forces acting on a suitable finite wing of the same airfoil section 
moving linearly with the same relative velocity and the same geometric 
angle of attack that is experienced by the propeller blade section. The 
contributions from all the blade elements may be added up using a strip 
integration technique to get the overall propeller characteristics. From 
the geometry of Figure 1 it can be shown that the incremental thrust and 
torque of any blade element operating at a geometric angle of attack, 

a = 3 - are given by 
S 

= \ BCpIV 2 (C L cos c 1 > - C D sin <f>) (2) 

and 

^ - j BC rpW 2 (C^ sin <f> + C D cos <j>) 


(3) 
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2 2 2 2 

where W = V + ft r . The validity of equations (2) and (3) is largely 
dependent on the characteristics of the wing used to represent the 
propeller blade section In order to accurately account for the inter- 
ference effects, the characteristics of the assumed wing would have to 
change from station to station However, the blade element theory greatly 
simplifies the problem by assuming the wings used to represent the blade 
element sections have an aspect ratio of six. Thus the problem is 
completely specified by equations (2) and (3) when the chord and blade 
angle distribution and the aerodynamic characteristics of an aspect ratio 
six wing of the same airfoil section as the blade element are known. 

The blade element theory better takes account the effect of the 
geometric shape of the propeller than does the momentum theory. However, 
it still fails to accurately account for the interference velocities 
induced by the trailing vortex system. 

i 

Vortex Theory - The vortex theory of propellers is basically a com- 
bination of the momentum theory and the blade element theory. Therefore, 
it is often referred to as the blade element-momentum theory; To avoid 
confusion with other theories it will be referred to as the Glauert vortex 
theory [5, 6, 7] m this paper. The Glauert theory is based on the 

assumption that the trailing vortex filaments which are produced by the 

* 

rotating blades form helical vortex sheets as they pass downstream. An 
exact application would require the induced velocity produced by the 
trailing vortex system to be computed at each blade station. However, 
Glauert simplified the problem by assuming the propeller to have an 
infinite number of blades. This assumption removes the periodicity of 
the flow, therefore, the velocity for any given radius is constant over 
the propeller disk. As a result, the momentum theory may be used to 
evaluate the interference velocities. The geometry of the problem is 
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given m Figure 2 a and a' are defined as the axial and rotational 

interference factors respectively V(1 + a) is the forward velocity, V, 

plus the axial interference velocity. flr(l - a') is the angular velocity 

of the propeller blade section minus the rotational interference velocity. 

The effective velocity, W , of the blade element is just the vector sum 

of V(1 + a) and f2r(l - a' ). From the geometry given m Figure 2 the 

following relations may be defined 

» 

W sm <j) = V(1 + a) (4) 

W cos <j> = fir(l - a* ) (5) 


It is often convenient to define a non dimensional speed ratio. A, 
as the forward speed divided by tip speed. Using equations (4) and (5) 
the speed ratio may be written m terms of the interference factors as, 

x ■ 5 r ■ TT-r tan ♦ (6) 

The equation for the elemental thrust and torque are found m the same 
manner as for the blade element theory, except now they are written m 
terms of the interference factors and are given m non-dimensional form as, 


R ■ = a ~ (1 - a') 2 (c cos <j> - C D sm 4>) sec 2 4> 

and 

d^c r ^ 2 2 

Cl - a') (C L sm 4> + C cos <j>) sec 4> 


(7) 

( 8 ) 


where the non dimensional thrust and torque coefficients are defined by 

T = T c tt R 4 p ft 2 (9) 

and 

Q = Q c tt R S p fi 2 (10) 

respectively Equations (6), (7), and (8) are not sufficient to determine 
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the interference factors a and a*. It is now necessary to make use of 
the assumption of an infinite number of blades . This makes it possible to 
apply the momentum theory to obtain relations needed to evaluate a and a' . 
The results of the axial momentum theory, where rotational momentum is 
neglected, can be used to determine the axial interference factor as 


o (C. cos <J> - C n sin 40 
a __ L D 

1 + a 4 sm 4> 

t 


(ID 


Equation (11) is not exact. In addition to neglecting the flow periodicity 
due to a finite number of blades, it also is based on the assumption of 
no slip stream contraction, and no rotational velocity component. The ro- 
tational interference factor may be evaluated by equating the torque to the 
rate of increase of angular momentum. This relation may be written in terms 
of the rotational interference factor as. 



tf(C sm <j> + Cy cos 40 
4 sm <f> cos <f> 


( 12 ) 


Equations (6), (7), (8), (11), and (12) are sufficient to determine 
the thrust and torque of a given propeller when the number of blades, 
the chord distribution, blade angle variation and the aerodynamic charac- 
teristics of the airfoil used at each blade location are known. 

Due to the assumptions made, the above set of equations accurately 
apply only to a lightly loaded propeller with a large number of blades. 
However, reasonable results may be obtained for many conditions varying 
greatly from the above restrictions The predicted thrust and torque 
distributions along the blade differ from the actual distributions for 
propellers operating at high blade loadings or for propellers with square 
tipped blades In the case of square tipped blades, the vortex theory 
fails to predict the fact that the thrust must fall to zero at the tip [8]. 
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Goldstein's Theory - Goldstein [9] dispenses with the assumption that 
the spacing between successive trailing vortex sheets is small and assumes 
there is a bound vortex filament springing from each blade element The 
local strength of the trailing sheet is the negative of the circulation 
gradient around the corresponding blade section The circulation must fall 
to zero at the root and the tip These trailing vortex filaments approxi- 
mately follow a helical path m the slip stream All the filaments together 
form a helicoidal surface. Goldstein specifically considers the optimum 
circulation distribution for which the energy lost due to the production of 
the trailing vortex system is a minimum, for a given thrust. Betz proved 
for lightly loaded propellers that this optimum circulation distribution 
corresponds to the requirement that the vortex sheet move rearward as if it 
were a rigid helicoid. Thus, the flow m the wake between the helicoidal 
surfaces is that of an mviscid, continuous, irrotational fluid, with zero 
circulation The circulation around a given blade element is equal to the 
discontinuity m the velocity potential of the respective point on the 
helicoidal vortex surface Since the problem is one of potential flow, it 
is possible to apply Laplace's equation, 

V 2 <J> = 0 (13) 

where <j> is the velocity potential. Using the relations 

C = 9 ~ ~ C 14 ) 

and 

r 


V 


V 


(15) 



Laplace's equation in polar coordinates may be transformed to a more 
convenient form given by. 
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U 


2 ,9* 


9y 


2 + ^ + + V) 


4 =0 


(16) 


The boundary condition for the problem becomes. 


8<j> _ _ wV 

9c " vWr 2 a 


( 17 ) 


for ? = 0 or jt and 0 <r <R. In addition, 4> must be continuous every- 
where except at the helicoidal vortex surface, and its derivative must 
vanish when r is infinite. Goldstein goes through a rigorous solution 
of this problem to obtain the loading factor, k, which is just a non- 
dimensional circulation defined by 


k = 


r b 

2w ir r tan cf) 


(18) 


Lock [10] applied Goldstein's solution to a method for predicting the 
characteristics of an arbitrary propeller. The development of this method 
parallels that of the Glauert vortex theory quite closely, except that 
Goldstein's solutions are used to determine the interference velocities 
at the blades The equations that are obtained for the interference 
factors are 


1 + a 


2, a C, cos <b - CL sin d> 
cos <j> L D 

k , 2 , 

4 sin <p 


(19) 


and 


. 2 , a C_ cos d) + C T sm (b 

a' _ cos d> D L 


1 - a' 


4 sin <f> cos <p 


( 20 ) 


bquations (19) and (20) differ from equations (7) and (8) of the Glauert 

2 . 

pQC (h 

vortex theory by the factor — ~ - It should be noted that m the limit 
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when the number of blades goes to infinity. 

k = cos^ cf> (21) 

Thus, equations (19) and (20) reduce to the Glauert vortex theory, 
equations (7) and (8), for the case of an infinite number of blades. 

Goldstein's loading factor, k, only applies exactly when there is 
a particular distribution of circulation along the blade to give optimum 
loading as defined m Goldstein's original problem. However, Lock shows 
that equations (19) and (20) make adjustments for tip losses and are 
generally more accurate than the Glauert vortex theory formulae, especially 
in the case of a square tipped blade. 

Theodorsen' s Theory - Theodorsen [11, 12, 13, 14, 15, 16] attacks the 
same problem of the optimum propeller as was developed by Goldstein, How- 
ever, Theodorsen expanded and extended the results of Goldstein's work 
Betz's requirement that the trailing vortex sheet move rearward as a rigid 
vortex sheet for an optimum loaded propeller was proved only for a lightly 
loaded propeller Theodorsen extended the interpretation of these results 
to apply to heavy loadings . 

Goldstein's original solution also applied only to infinitely light 
loadings Theodorsen showed that it is possible to extend these results 
to heavy loadings by including the rearward displacement velocity of the 
helicoidal vortex surface and doing all the calculations in the wake 
infinitely far behind the propeller This does not change the potential 
problem that must be solved. However, it does change the operating con- 
ditions which a given solution corresponds to. Including the rearward 
displacement velocity of the vortex surface, instead of assuming it to be 
negligibly small, takes account of the contraction of the slip stream. 



Theodorsen's method depends on the solution for the propeller loading 
coefficient, k, which is a function of the distribution of circulation 
along the blade. Goldstein's results may be used to evaluate the loading 
coefficient. However, Goldstein was only able to obtain a small number 

n 

of solutions. Iherefore, Iheodorsen used an electrical analogy method 
to experimentally measure the loading function By using the electrical 
analogy, he was able to obtain the loading function for dual rotating 
propellers as well as for single propellers. The drawback of this method 
is that one is limited to the charts and tables which Theodorsen gives 
for the loading function. 

Application of the Prandtl Theory - Lock [17, 18, 19, 20] put forth 
a theory which solved the propeller problem through an application of the 
Prandtl theory [21] for a monoplane wing In calculating the induced 
velocity along the blade, the actual helicoidal trailing vortex sheets are 
replaced with a system of vortices m two-dimensions. This method 
assumes that the flow around the helicoidal vortex sheets can be approxi- 
mated by the flow around a system of parallel two-dimensional vortex 
sheets whose strength is equal to the actual strength of the helical 
vortices at the same radial distance from the axis of rotation. The 
distance between successive sheets is equal to the pitch of the airscrew 
The Fourier series may then be used to represent any arbitrary distribution 
of vorticity with radius Strict application of this theory would require 
the coefficients of the Fourier series to be found in a manner analogous to 
that employed in wing theory, where a system of simultaneous equations m 
the coefficients must be solved. However, to simplify the technique, a 
graphical method is applied to determine the coefficients of the series and 
corrections are then applied to account for errors introduced by replacing 



the helicoid vortex sheet by a system of two-dimensional plane laminae 
The method m its final form is strictly a graphical technique, using charts 
developed by Lock 

Morxya Theory - Moriya [22] developed a calculation method for the 
aerodynamic characteristics of a propeller m which he obtained the induced 
velocity by introducing an induction factor, I To obtain the induction 
factor it is assumed that the trailing vortices springing from each blade 
form a helicoidal surface which extends rearward from the blade to infinity 
The downwash velocity at each blade element due to the whole vortex system 
is then calculated by application of the Biot-Savart law The downwash 
integral becomes singular when calculating the velocity induced at a point 
where the trailing vortex filament springs from that same point. If 
the vortex filament were straight instead of helical, it would exhibit the 

same singular behavior Therefore, Moriya introduced an induction factor, 

I, which defined as. 


I = 


d w 

n 

d w 


n l 


( 22 ) 


where d w is the normal velocity induced by a helical vortex filament and 

d w is the normal velocity induced from a straight vortex filament At 
n l 

the point where the vortex filament leaves the blade, the induction 
factor takes on the value of unity Induction factors for various advance 
ratios were calculated graphically by Moriya and are tabulated m reference 
[ 22 ] 

A successive approximation scheme to determine the aerodynamic char- 
acteristics of an arbitrary propeller was developed by Moriya In this 
scheme the circulation around each blade element is calculated corresponding 
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to an equivalent monoplane wing of aspect ratio six Using this calculated 
circulation distribution and the suitable induction factors, the normal 
velocity at each station is obtained by graphically integrating the following. 


w = 
n 



(23) 


An equivalent aspect ratio for each blade section can then be calculated 
from 




(24) 


This scheme is repeated by calculating the circulation distribution 
corresponding to the new equivalent aspect ratio at each blade element 
This iteration is repeated until the equivalent aspect ratios of two 
successive approximations are sufficiently close The lift and drag coef- 
ficients are found m the final iteration and are used m the normal strip 
theory formulas, equations (2) and (3), to find the thrust and torque 
coefficients It is reported by Moriya that only one or two iterations are 
needed to obtain a good degree of accuracy 
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III THEORETICAL INVESTIGATION 


Most of the preceding theories involve many simplifying assumptions 
to obtain the interference velocities needed to analyze the aerodynamic 
characteristics of a propeller The equations and method needed to solve 
for the induced velocities with a minimum of simplifying assumptions m 
a manner analogous to that of finite wing theory will now be developed 
. The basic assumption is that of a helicoidal trailing vortex sheet which 
extends rearward from the blade to infinity. The velocity induced by a 
segment, d £, of a trailing vortex filament can be calculated at each blade 
section by means of the Biot-Savart law which is given as. 


5 v (r) = 


T d& x (t-s) 
4 tt 


(25) 


|t~s | 3 

where t is the vector from the origin to the blade section at which the 
induced velocity is to be calculated and s is the vector from the origin 
to the' vortex filament segment dZ. For the geometry depicted m Figure 3 
the following relations can be wnt’ten. 


t = R x^i . (26) 

/ A A A 

s’ = R x cos 0 i + Rx sm 0 j + Ra tan <j> k (27) 

A A A 

d£ = R x sec <j> d0 [- cos <p smdi + cos <j> cos 63 + sm <j> k] (28) 

The induced velocity at a point x^due to a semi-infinite vortex filament 
springing from the point x on the blade is found by substituting equations 
(26), (27) and (28) into equation (25) and then integrating this equation 
over 0 from zero to infinity The normal component of this induced velocity 


is given as. 


dw = 
n 


T 

4irR 


/. 


“ Ir 2 

x| X 


xx^~ cos 0 ] + [x(9 sm 0 + cos 8 ) - x'*] 


[x 


_2 


2 xx ■* cos 0 + 


X 2 + 


0 ' 


x 2]3/2 ^2 + 'Zjl/Z 


d 0 


(29) 
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The normal induced velocity at the point x' due to the vortex filaments 
springing from the radius x on all the blades is 


df 


B 


* 

[x - xx' cos 9^]x / + [X(0 sm 0^ + cos 0^)-x"] A d 


dW " ' Si I [x' 2 - 2XX' cos 6 ♦ x 2 X 2 ] 3 ' 2 [A 2 T X'] 1/2 


d0 (30) 


i a n ^ 2TT (B-k) 

where 0. = 9 + £ — 

k B 


The normal induced velocity at some point x" due to all the vortex filaments 
springing from each of the blades is found by integrating equation (30) over 
the blade and can be expressed as 


» n (x'J - 


<ir . b 

dx 


4irR 


*E f ^ 
k=1 Jo [X 


^-xx" cos S^Jx^ [x(0 sm 0j, + cos 0^)-x'" A^ 


- 2xx'" cos0. 


2 A q 2-, 2,3/2 2 , 2,1/2 

+ x +0X] [X +x^j 


d0dx 


(31) 


Equation (31) is quite complex and is difficult to handle since the inner- 
most integral contains a singularity at x = x' that is not shown explicitly. 
This singularity can be made explicit by applying the induction factor, I, 
introduced by Moriya [22] and defined by equation (22) The induction 
factor is simply the ratio of the normal induced velocity for a helicoidal 
trailing vortex system to the normal induced velocity for a straight trailing 
vortex system of the same strength The induced velocity has a component 
from the vortex sheet springing from each blade The induced velocity for 
the helicoidal vortex sheet springing from the blade at which the calcula- 
tions are being made becomes infinite at x = x' The induced velocity for 
a straight vortex system also becomes infinite at x = x" The infinite 
, velocity m both cases is a result of that portion of the vortex filament 
m the neighborhood of the point x' . Since the property of the infinity at 
x = x' is the same in both cases, the ratio of the two induced velocities 
must be unity at the point x = x". The contributions to the induction 
factor at x = x' from all other helicoidal trailing vortices springing from 
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all other blades must be zero since only the blade at which the induced 
velocity is being calculated has the property of infinity at x = x'*. There- 
fore the induction factor takes on the value of unity at x = x^. 

The induction factor, which is a continuous function of A, x and x^ 

i 

only, may be evaluated numerically by computer. When x approaches x" the 
induction factor becomes hard to numerically determine with accuracy. 
However, since the induction factor is a continuous function and equals 
unity when x = x' , the values when x approaches x' may be interpolated 
to a high degree of accuracy Typical curves of the induction factor 
for A = 0.3 are shown m Figure 4 

The normal induced velocity at a point x' may now be expressed 
compactly as. 


v ( x ") = 


1 dr 

dx I 


-x dx 


( 32 ) 


4ttR x-x 

The singularity at x = x" is now shown explicitly m equation (32) 

The method is now formulated by applying the fundamental equation of 
finite wing theory. 


a = a - a 
g o I 

a is the geometric angle of attack and is defined by 
§ 

a = 3 - Arctan (^) 
g x 


(33) 


(34) 


The effective, or sectional angle of attack, a Q , is evaluated by applying 


the Kutta-Joukowski theorem and is evaluated as. 


2 r 


R SI (A 2 + x" 2 ) 1/2 a -C 
k J o 


The induced angle of attack, a is given by 


w w 

„ _ n • n 

a i = tan IT = W 


(35) 


(36) 




A- 23 

Putting equation (34), (35), and (36) into equation (33) results m the 
fundamental integral equation which must be solved for the circulation, T 
The boundary condition for the circulation requires that, 

T(0) = T(l) = 0 (37) 

Therefore, we may express the circulation as a sine series given by 

CO 

rex') = E A sin mix'* (38) 

n=l 

Using this series, the fundamental integral equation can be expressed as 


2EA sin n ttx' 
n 


l n irA 


n 


a c R tl (A 2 + x' 2 ) 112 
o ^ 


4tt R 2 ft (A 2 + x' 2 ) 1/2 


/A 

COS n _TTX 


X - X 


g - Arctan (— *•) 
K x 1 


(39) 


Difficulty now arises m determining the second term m equation (39) since 
it contains a singularity m the integrand at x = x' This problem is 
handled by breaking the integral into three separate integrals, 


J l 6 ~x'+& 

M dx = I M dx + I M dx + I M dx 

o •'o J x'-6 'x'-S 


(40) 


where M = 


o 

cos nirx 


I and 6 is a small but finite value The first and 


x - x 

third integrals on the right-hand side* of equation (40) may be evaluated 

using established numerical techniques The second term, however, requires 

special consideration It can be shown that the integrand, M, will go to zero 

-x+<5 


m the limit as 6 goes to zero Thus, 


/•A * v. 

. x-d 


M dx will have a finite value. 


A method to evaluate this integral for a finite 6 is required This is 
accomplished by first making a variable substitution using, 

- x - x" (41) 


T 
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The integral can then be written as 

,6 


L 


COS n Tt (t + x') 


I dT 


(42) 


Equation (42) is then expanded m a power series m t It should be noted 
that the induction factor, I, is a summation of integrals. The induction 
factor is a continuous function by virtue of the manner in which it was de- 
fined However, each of the integrals making up the induction factor is 
not continuous and these integrals do m fact contain smgulantes There- 
fore, it is necessary to expand the induction factor in a power series m a 
manner which will eliminate these singularities. Since the singularities 
occur for the helix angle, 0, equal to zero, they may be eliminated by in- 
tegrating the singular integrals from 0 equal e to 00 and then separately ex- 
panding the integral from 0 equal 0 to e, which contains the singularity. 

Once the power series m T is obtained, it is integrated term by term and 
2 

terms of order t and larger are dropped This gives the following approxi- 


mate relation. 


x+6 


cos n ttx 


x-6 


I dx = 25 cos (nfix^) ^ 

k=l J 


2 f x '+ 6 

M dx + 26 


X"~ 5 


x"cos (mrx") 


on 2 -2. 
2(A +x ) 


(1 + £n(2e 2 (A 2 + x"* 2 )) - nff 2 sm mrx' 


6£n (6) 


dx^cos (niTx'') 


? 2 
2 (x + A j 


(43) 


where e is a small value and e » 6. The singular integral can thus be 
approximated for the finite region of width 6 around the singularity. 

In order to evaluate the coefficients, A , of the infinite sine series 
given m Equation (38), the series is truncated at N terms. The integral 
Equation (39) is then applied at N locations along the propeller blade. 

This gives a system of N equations m the N unknowns This system of 
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equations is solved for the coefficients, A^, which are used to evaluate 
the circulation, using Equation (38) The effective angle of attack can 
then be found using Equation (35) . Once the effective angle of attack is 
found at each designated position along the blade, the lift and drag co- 
efficients can be found directly from the two-dimensional aerodynamic 
characteristics of the airfoil used at that position on the propeller 
blade The differential thrust and torque coefficients can be found by 
applying Equations (2) and (3) . 

The preceding method has been programmed for use on a digital com- 
puter. An explanation of the program is given m Appendix A and a list- 
ing of the program is given m Appendix B. 



A-26 


APPENDIX A 

Computer Program Explanation 

The computer program listed m Appendix B was developed following 
the method proposed m Section III of this report The variables in 
Table V must be inputs to the program. 

Table V Input Variables to the Computer Program 
AO - Lift curve slope per radian. 

BT - Blade pitch angle m degrees which must be specified at N locations 
along the blade. 

C - Non-dimensional chord which must be specified at N locations along 
the blade. 

CDO, CD1, CD2 - Drag coefficients of the airfoil section where the total 
drag coefficient, CD = CDO + CD1 a + CD2 a^. 

DELI - Interval around the singularity which is not integrated directly. 

EL - Tip speed ratio. 

N - Number of terms m the sine series 
NB - Number of blades. 

R - Radius of the propeller in feet. 

SKPIl, SKPI2 - The upper and lower limits, respectively, around each of 
the N blade locations, between which the induction factor 
is interpolated 

WO - Rotational component of tip speed of the propeller in feet per second. 
X,W - Languerre- Gauss integration coefficients. 

ZTP - N locations along the blade at which the calculations., for the sine 


series are made 
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The program uses a 44 point Languerre-Gauss integration scheme [25] to do 

the integration from zero to infinity in order to evaluate the induction 

* 

factors. Induction factors near x’ are Interpolated using a spline Inter- 
polation formula [26] 



where the C's are the spline coefficients generated by the spline inter- 
polation routine developed in reference [27] . The spline integration 
technique has been found to be very accurate since it fits a curve with 
continuous derivatives through second order to the given function and then 

integrates the area under this curve. This integration scheme is also used 

* 

in calculating the total thrust coefficients from the differential thrust 
coefficients . 
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to become less accurate due to the integrand M becoming very large. A 
value of 6 between 10~ s and 10~ G has been found to give good results. 

The matrix equation for the coefficients of the sine series is solved 
using Gauss elimination. The answer is checked for accuracy by substituting 
it back into the matrix equation. 



€3'' „ 

f=-»] 

gg 

Cf 

£> P 5 * 

F S 

hj O 


Go 


n 


C THIS PROGRAM COMPUTES THE THRUST COEFFICIENT FUR AN ARBITRARY 

C PROPELLER DESIGN USING THE METHOD DEVELOPED IN CHAPTER THREfc 

C OF THIS PAPER. 

C 

C 1 

C IMPUT VARIABLES TO THE PROGRAM ARE: 

C 

C AO - THE LIFT CURVE SLOPE PER RADIAN 

C BT - THE BLADE ANGLE AT EACH OF THE N LOCATIONS ALONG THE BLADE 

C CD 0 r GDI, CD2 - SECTIONAL DRAG COEFFICIENTS 

C C - THE CHORD NON-O I MENS I ONALI ZED ^ I TH THE RADIUS 

C DELI - THE INTERVAL AROUND THE SINGULARITY THaT IS INTEGRATED 

C SEPARATELY 

C EL - THE TIP SPEED RATIO = V/WG*R 

C EP - A SMALL NUMBER WITH EP»DEL1 

C N - THE NUMBER QF POINTS IN THE SINE S-ERIES 

C N8 - THE NUMBER OF PROPELLER 8LADES 

C R - THE RADIUS OF THE PROPELLER IN FEET 

C SKPI1»SKPI2 - BLADEWISE LOCATIONS BETWEEN WHICH THE INDUCTION 

C FACTOR 'IS INTERPOLATED 

C WO - THb ROTATIONAL SPEED IN RADIANS/SEC 

C X»W - LANUUERRE-GAUSS INTEGRATION POINTS 

C ZTP - THE N NON-DIMENSIONAL BLADEWISE LOCATIONS AT WHICH 

C THc CALCULATIONS ARE MADE FOR THE SINE SERIES 

C 

IMPLICIT REAL*8 (A-H.O-Z) 

INTEGER-* I0UMI9 i 

DIMENSION ZTPI81), ZTI81), XZ(3L), Q181), AA ( * » 8 1 ) » SSI 5* 81 ) , 

1 GINTI9 ) , S I NT ( 9 ),Slol), RINTI9 ),ALP019 ),C(9 ), 

2 YI(81), CA(9 ,9 ), RH S ( 9 ), ASI9 I, CC A ( 9 ,9 ), 

3 I S ( 3 1 ) j TIGRU IL , 11) , DTU1L), BETA! 11) 

COMMUN/BBB/XI 44) , W(44) 

COMMON /CCC/ TT t TTP , PI, EL, EP, NB 
C RADIANS TO DEGREES CONVERSION 

RD„C = 57. 295779513082300 
PI=3.14159265358979D0 

C READ IN DATA FOR A *4 POINT L ANGUERRE-GAUS S INTEGRATION 

RE AD ( 5 , 5 ) t X ( L ) , WIL), L=l,44) 

5 FORMAT (2D30.0) 
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WRITE (6,6) ( X ( L ) , W ( L ) , L=l,44) 

6 FORMAT (« », 030*15, 030.15) 

N= 9 

NB = 3 

EL=O.2906DO 
EP=0. 3000 
W0=153. 93300 
A0=O.07D0-RDC 
R= 1 . 4000 

CDO=O.0333432lDO 

CD I = — 0 .00800893 DO *ROC 

CD 2=0. 0006339300^00 *=* 2 

C NP IS THE NUMBER OF POINTS USED IN THE SPLINt INTEGRATION 

C NP MUST BE AT LEAST 4*N SO THAT WHEN INTEGRATING FRON X=0 TO 

C ZTP-DEL OR ZTP+DEL TO 1 THERE WILL ALWAYS BE AT LEAST 4 POINTS 

C FOR THE SPLINE INTEGRATION. (SUBROUTINE SLP1Z WILL FAIL WITH 

C LESS THAN FOUR POINTS/ 

NP={N-L)£10+1 
DO 9 J = 1 ,N 

C READ IN, THE BLAOE ANGLE, CHORD, AND CORRESPONDING LOCATIONS 

C ON THE BLADE 

RE AD (5, 7) BT, C(J), ZTP(J) 

7 FORMAT (3F10.6) 

8ETA( J )=BT/RDC 

C AL PC IS THE GEOMETRIC ANGLE OF ATTACK 

IF ( ZTP ( J ) . L E.O .000 LDO ) AL PO ( J ) = B T/ RD C-P I /2 .0 DO 

IF ( ZT P ( J ) . GJ .0 .QOO IDO ) ALPO(J)= B T/ ROC-DA T AN l EL / Z TP ( J ) ) 

WRITE 16, d) ZTP(J), BT, ALPO(J), C ( J ) 

8 FORMAT ( » * , 4F20.6) 

9 CONTINUE 

DO 10 1=1, NP 

10 ZT ( I ) = FLOAT ( I-IJ/FLOAT (NP-1) 

OEL1=O.OOOOG1DO 

DO 400 J=1 ,N 
DO 4 1=1 , N P 
IS I I ) = 0 
4 Y I ( I ) = 0. ODO 

WRITE (6, LI) ZT P { J ) 

11 FORMAT (*0*, *ZTP=», F8.3) 
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SH=EL**2+ZTP ( J)**2 


THIS SECTION CALCULATES ALL THE INDUCTION FACTORS NEEDED 
TTP=ZTP{ J) 

tt=ttp 

NOPT = 2 

PI2 IS THE SUM OF THE INDUCTION FACTORS FOR BLADES I AND 2 
DIVIDED BY ( ZT-ZTP ) 

CALL CALCI { P 1 2 , NOPT) 

NOPT = 3 

P3 IS THE INTEGRAL FROM EP TO, INFINITY OF THE INTEGRAND M 
WHERE I3=( ZT-ZTP ) ^INTEGRAL OF M 
CALL CALCI ( P 3 t NOPT) 

NOPT=l 
TT P= ZTP I J ) 

NAB = 0 
NCT K=0 * 

ICTR=0 

C INTERPOL-ATE VALUES OF I BETWEEN SKPI1 £ SKPI2 

C WHERE SKPIL £ SKPI2 ARE READ IN FOR EACH ZTP LOCATION 

RE AD ( 5 » 1 2 ) SKPI1, SKPI2 
L 2 FORMAT (2FL0.3) 

WRITE (6*13) SKPIi, SKPI2 

13 FORMAT MO*,' SKIP DIRECT I CALCULATION BETWEEN ZT=*, F7.3, 

1 * AND' , F7.3 ) 

DO 15 I = 1 * N P 
TT = ZT ( I ) 

c FIND THE VALUE OF THE DO LOOP PARAMETER FOR WHICH ZT=ZTP 

IF (DABS( TTP-TT) . LT . 0 . 00 D 00 1 DO ) ICTR=I 
IF {TT.LT. SKPII. OR. TT.GT.SK D I2) CALL CALCKXI, NOPT) 

C STORE THE INDUCTION FACTOR IN THE ARRAY YI 

IF ( TT.LT. SKPII. OR. TT.GT.SKP 12 ) Y I ( I )=»X I 
C STORE WHICH INDUCTION FACTOR CALCULATIONS HAVE 8EEN SKIPPED 

IF (TT.GE. SKPIL. AND. TT.LE.SKPI2) ISII) = 1 
15 CONTINUE 

DO 17 1=1, NP 
ICK=0 
TT=ZT ( I ) 
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o o n n o on 




IF ( IS II ). NE. 1. AND. TT.GE.SK PI 1-0. 074D0. AND- TT. L£ . SK P I 2 +0 . 074D0 
1 .OR. I .EQ. ICTR) ICK=2 

NAB IS THE NUMBER OF POINTS AT WHICH THE DIRECT CALCULATION 
OF THE INDUCTION FACTOR HAS BEEN CALCULATED 
IF (ICK.EQ.2) NAB=NA8+1 
IF (ICK.EQ.2) XZ(NAB)=TT 

STORE THE CGMPUTEO INDUCTION FACTORS IN ARRAY Q 
IF ( ICK.EQ.2. AND. I. NE- ICTR ) Q(NA6)^YI(I) 

IF ( I.EQ.ICTR) Q( NAB) =1.000 
17 IF (I.EQ.ICTR) NCTR=NA8 

WRITE ( t> » 2 1 ) NAB, ICTR, NCTR 
21 FORMAT ('O', * NAB- * , 3110) 

WRITE (6,19) ( XI l L ) , Q(L), L= 1 , NAB ) 

19 FORMAT ( 1 • , 2F20 .5 ) 

USE THE SPLINE INTERPOLATION FORMULA TO OBTAIN THE INDUCTION 
FACTORS THAT WERE SKIPPED IN THE DIRECT CALCULATIONS 
SLP1Z IS A COMPUTER SUPPLIEO SUBROUTINE WHICH WHICH GENERATES 
SPLINE INTERPOLATION COEFFICIENTS. 

CALL SPL1Z (XZ, Q, NAB, AA, SS ) 

DO 18 1=4, NP 
TT = ZT ( I ) 

C USE SPLINE INTERPOLATION FORMULA 

IF ( IS ( I I .EQ. 1 . AND. I ,LT. ICTR) 

1 Y I ( I } = A A ( 1 ,NCTR~l)*(XZ( NCTR )-TT 

1 +AA(2 ,NCTR-l)* (TT-XZ(NCTR-i) >*-*3 

2 +AA( 3,NCTR-1)* (XZ(NCTR)-TT) 

3 +AA(4,NCTR-1)«(TT-XZ(NCTR-1) ) 

IF { I .EQ. ICTR) Y I ( I ) - 1 . ODO 

IB IF ( IS ( I 1 . EQ. I . AND. I .GT. ICTR ) 

L Y I ( I ) = AA ( 1 ,NCTR )*( XZ (NCTR + U-TT ) 

2 +AA(2,NCTR)*-(TT-XZ(NCTRn#*3 

3 FAA( 3,NCTR ) * < XZINCTR+D-TT) 

'4 +AA(4,NCTR)v(TT-XZ(NCTR) ) 

IF I ZTP( J) -GT.DfcLL) TT=ZT P l J ) -D L LI 
IF (ZTP( J) .GT.DELL) 

1 XI BZ = AA( I , NCT R- 1 )*(XZ(NCTR)-TT)**3 

1 +AA( 2 , NCTR- 1)4 (TT-XZt NCTR— I ) )**3 

2 +AA(3,NCTR-L)*(XZ(NCTR)-TT) 

3 +A A ( 4 , NCTR— 1 ) * { TT— XZ I NCTR— 1 ) ) 
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IF (1.000- ZTP (J ) . GT . DEL L ) TT = ZT P ( J ) +DEL 1 
IF ( l.OOO-ZTP(J) .GT.DEL1 ) 

1 X I AZ= AA ( 1 , NC TR )*(XZ(NCT R+ 1 ) -TT ) 3 
Z +AA(2,NCTR)*(TT-XZ(NCTR) )*-*3 

3 + AA(3,NCTR)*(XZ(NCTR+-L) — TT) 

4 + AA ( 4 , NCTR ) *( TT-XZ ( NCTR) ) 

C ALL THE INDUCTION FACTORS HAVE NO* bEEN CALCULATED AND 

C STORED IN VI, XI6Z, XIAZ. XI3Z & X I AZ ARE THE INDUCTION 

C FACTORS AT ZTP-DEL AND ZTP+OEL RESPECTIVELY. 

D'J 210 NS = 1,N 
APG=NS*P 1^ ZTP ( J ) 

WRITE ( 6 , 20 > NS 

20 FORMAT (’O', *+++ NS = ‘, 13, 1 4- + -I- + + ++ + + + +• ) 

DO 22 I J = 1 , N 
QINT{ I J ) =0 .0 DO 
RINT( I J ) — 0. ODO 
22 S INT ( I J )=0.0D0 
DE L = OE L 1 
N0 = 0 

IF (ZTPf-J).LE.DEL) GO TO 53 
C 

C THIS SECTION CALULATES THE INTEGRAL FROM 0 TO ZTP-DEL 

C 

DO 25 1=1, NP 
ND=ND+ 1 

IF (ZT(l )-Z T P(J) .GE.-OELI ) GO TO 30 
XZ ( I > = ZT { I ) 

C STORE THE INTEGRAND IN ARRAY 0 

25 Q( I ) = OCOSlNS*PI*ZT< I ) )*YI ( I) / ( Z T (II -ZTP ( J ) ) 

30 XZ{NO)=ZTP(J)-DEL 

U(ND} = DCOS(NS- i< PI*XZ(ND ) )*XI3Z/ (XZ(ND)-ZTP(J) ) 

ND1=ND-L 

TF (NS.EQ.l) WRITE (6,35) (XZ(L), Q( L ) , YI(L), L=1,ND1> 

IF (NS.EQ.l) WRITE (6,35) XZ(NU), Q ( ND ) , XI6Z 

35 FORMAT ( ' ’ , 3D15.5 ) 

C INTEGRATE Q USING SPLINE INTEGRATION FORMULA 

CALL SPL1ZIXZ, Q, ND , AA, SS ) 

SUMG=0.0D0 
ND 1 = ND— 1 
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DO 40 K=1,ND1 

40 SUMQ=SUMQ + .25D0*AA{l,'k}*(XZU<+U-XZ(K) ) **4 

1 + .2 5D0*AAl 2 ,K)*(XZ( K+i l-XZ(K) )**4 

2 + .50D0*AA(3,K)*(XZl,K*-l)-XZ(i<) )**2 

3 +.50D0*AA(4,K)*(XZ(K+L}-XZ{K) )**2 
Q I NT (NS) -SUMO 

WRITE {6,50) SUMQ 

50 FORMAT {’ ', 'INTEGRAL OF Q= ' , Eli. 3) 

53 CONTINUE 

IF ( IGT.Ed.l ) GO TO 75 

IF (ZTP(J).LT.DEl.OR.UODO— ZTP(J).LT.DEL) DEL=DEL/ 2.0D0 
THIS SECTION CALULATES THE INTEGRAL FROM ZTP-DEL TO ZTP+DEL 
ATERM=DC OS ( ARG ) *2 . 0D0*DEL*Pi2 

BTERM1=2.0D0*DEL*( (ZTP ( J )*DCOS( ARG) ) / l DSwRT ( 2- ODD ) *SH ) 

1 *<l.ODO+DLOG(2.0DO*EP**2*SH) > 

2 -NS*Pl*DSQRT{Z.ODO)*DSIN(ARG) ) 

BTERM2* -DEL* OLOG( DEL L }*( 4-0D0*ZTP < J ) *DCOS ( ARG) 

1 * /( DSQRT (2 .0D0)*SH> ) 

CTERM=2.0D0*DCOS( ARG5*DEL*P3 
RINT(NS)=ATERM+BT ERM L +-BT ER M2 +CT EkM 
IF (NS.EG.l) WRITE (6,55) PL2, P3 
55 FORMAT {*0', 2E20.5) 

IF (NS.EQ.l) WRITE (6,60) A TERM, BTERM1, 3TERM2, CTERM 
60 FORMAT ( '0« , 4E20-5) 

WRITE (6,70) R INT (NS ) 

70 FORMAT (• ', 'INTEGRAL OF R=», Ell. 3) 

DEL=DEL1 

75 RINT(NS)=0.000 

IF < 1 . ODO-ZTP ( J ) . LE.DEL) GO TO 160 

THIS SECTION CALULATES THE INTEGRAL FROM ZTP+DEL TO 1.0 
ND=1 

XZ(ND)-ZTP ( J ) +DEL 

S(ND)=DC05{NS*PI*XZ(ND) )*XIAZ/(XZ(ND)~ZTP(J) ) 

DO 100 1=1, NP 

IF ( ZT { I l-ZTP ( J ) .LE.DEH ) GO TO 100 
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ND = N D + 1 
XZ (NO) =ZT( I ) 

C STORE THE I NEGR AND IN ARRAY S 

S(NO)=OCUS(NS*PZ*ZT< 1) l *Y1 < I I/(ZT( I)-ZTP( J) ) 

100 CONTINUE 
110 CONTINUE 

IF (NS.6Q.1) WRITE (c,35) XZ ( 1) , S(L), XIAZ 

LND=NP-ND 
DO 115 L=2 » ND 
LL = L + L NO 

IF (NS.EQ.1J WRITE I6,3:>) XZ(L), S(L), Y I ILL ) 

115 CONTINUE 

C INTGRATE S USING SPLINE INTEGRATION FORMULA 

CALL SPL 1 Z ( XZ * S, NO, AA , SS ) 

SUN S=0 . ODO 
ND L = ND-1 
DO 140 K = 1 , N D 1 

140 SUM S= SUM S+. 2 5D0*AA ( 1 , K ) * ( X Z ( K+ L ) -X Z ( K) )**4 

1 +.25DO*AA{2,K5*(XZ(K+l)-XZ(K) ) **4 

2 * + .50D0*AAI3,K»*(XZ(K+U-XZ(K) )**2 

3 +.50C0*AA( 4,K)*( XZ ( K+ 1 ) -X Z ( K ) ) **2 
SI NT { NS ) =SUMS 

WRITE l 6, 150) SUMS 

150 FORMAT ( ' 'INTEGRAL OF S= * , E12.4) 

C TINT IS THE TOTAL INTEGRAL FROM ZERO TO 1 

160 TIN T=0 INTINS) +R I NT (NS) + SINT INS ) 

T I GRL ( J , NS ) =T I NT 
WRITE! o, 170) TINT 

170 FORMAT ('O', 'TOTAL INTEGRAL FROM 0 TO I NF IN I T Y= ' , E20.5) 
C SET UP THE MATRIX EQUATION 

CA(J,NS)=2.0D0*DSIN(NS^PI !!: ZTP{J) ) / AO-NS^Cf J ) *T I NT /(h -.ODO-? 
CCA( J,NS )=CAt J ,NS) 

210 CONTINUE 

RHS( J ) = WO*R§§g£*OSQRT( SH}*ALPO(J)*C(J) 

WRITE (6,350) 

350 FORMAT ('0* , • 

400 CONTINUE 

WR I TE ( 6 , 500 ) ( ( CA ( M , L ) , L=1,N), M=1,N) 


se^v 



500 FORMAT (• » f 9FH.5) 

WRITE I 6 » 5 50 > IRHS(L), L=1.N) 

550 FORMAT ('O', 9F11.3) 

- SOLVE THE Matrix EQUATION FOR THE AS i/aiiicc lie IMP A 

CAL^GAUSZ^rA 1 ED GAUSS El ^INATION TECHNI QUE. 

ud\tc G ^ SZ ' CA ' AHS ' I OUM t IER) 

WRITE ( & , o 00 } I AS { L ) , L=1,N) 

600 FORMAT * r • a = ' , D20.5) 

WRI T £ (6 ,o25 ) 

625 FORMAT CO', UX, ! ’ 1 6X ’ 'GAMA', 16X, 'ALPHA I', 

DO 700 NK-l.N “ ° ' 14X ’ ' DTC '’ UX ’ -PHI ■ ) 

ST IGR=0. ODO 
SGAM=0 .000 

VC»hO*Rg^OSQRT{EL**2 + 2TP(NK)**2I 

VR IS THE RELATIVE VELOCITY 
V*=VC/ R 
DO 650 K=1,N 

ST I GR= ST IGR + T IGRL (MK , K ) * AS ( K }*K 

CALUL ATC THE CIRCULATION 

o5 0 SGAM»SCAM+AS(K)*DSIN{K*PI*ZTP{NK) ) 

AL I=STIGRM4.0D0*VR) 

CALCULATE THE EFFECTIVE ANGLE OF ATTACK 

ALG=2.0D0*SGAM/{ACJ*C(NK)vvC) 

PHI=BETA(NK)-ALO 

CALcStA^^^E^I^’c^fc 1 ;^^^ 410 

«L^^;r T -^ i o o s^sr^;c o ^r ROc,,ALo -°- io8ooi -- rou 

CD = CDO+CO 1 ■•’'ALO* COZ^A L0^*2 
CY = C L* DC OS ( PH [ )-CO*DS IN( PH I ) 

CALCULATE THE DIFFERENTIAL THRUST 

WRITE (6,6751 C ZTR pam /{ a * ?°°* P 1 2 * w0 ** 2 5 * ( 0 - 2 5 DO*P 1**3 ) 

675 FORMAT [- 6E20^) ’ AL 1 ’ AL °’ DTC ‘^>’ PHI 

700 CONTINUE 

FORMULA^ THC THRUSr 0VER THE 6LADE USING SPLINE INTEGRATION 
CALL SPLUUTP, DTC, M, AA , SSI 


)£-V 



o n o n n o o 


NM1=N- L 
SUMT=O.ODQ 
DO 702. K= l , NM L 

702 SUMT=SUMH-.250D0*AA{2,K)*(ZTP(K4-L)-ZTP(K))* ; i‘4 

1 -*-.5OGD0*AA{4,K) *( ZT P ( K + i ) -ZTP ( K ) ) # * 2 

2 + .2 50D0*AA{ ltlOM ZTP( K + LJ-ZTP (K) ) 

3 500D0*AAI3,K)*( ZTP(KM) -ZTP(K)) **2 
WRITE (6,705) SUM! 

705 FORMAT ('O', ‘THRUST COEFFICIENT^', F20.7) 

CHECK TU SEE IF THE MATRIX EQUATION WAS SOLVED CORRECTLY 
DO 725 1=1, N 
CKRHS=0-0 
DO 710 J =1 , N 

710 CKRH5=CCA( I , J ) *AS ( J) +CKRHS 

725 WRITE (6,730) CKRHS 

730 FORMAT ('O', F20.3I 
RETURN 
END 


vf 


SUBROUTINE CALCI (XI, NOPT ) 

THIS SUBROUTINE CALCULATES THE INDUCTION FACTORS AS FOLLOWS: 
NOPT = 1 CALCULATES ri+I2+I3=I 

NOPT =2 CALULATES THE INTEGRAL OF P1+P2 FROM 0 TO INFINITY 

mScIc = , 3 , G ALULATES THE INTEGRAL OF P3 FROM EPSILON TO INFINITY 
WHERE 1= (ZT-ZTP)*(Pi+P2+P3 ) 

IMPLICIT REALMS IA-H,Q-Z) 

COMMON / AA A/ KK, NNGPT 

COMMON /C CC / TT, TTP , PI, EL, EP, NB 

DIMENSION AN (3) 

EXTERNAL AUX 
NNOPT=NCPT 


DO 10 KK=I ,3 

IF (NOPT.EQ. 2- AND.KK.EQ.3) GO TO 20 
IF (NOPT .EU-3-AND.KK.NE. 3 ) GO TO 10 
CALL GLQU (AUX, ANS ) 

10 AN ( KK ) -ANS 

IF (NOPT -EQ.I. AND. DABS (TT-TTP) -LE. 0.00 IDO) 
i X I = I . ODO 


Z£--V 



ooooo noo 


i 


IF ( NOPT.EQ. 1.AN0.0ABS/(TT-TTP) .GT. 0.00 IDO) 

I X 1= l AN {(&) +AN ( 2 ) + AN (3 )) #( TT-TTP ) 

20 IF (NOPT.EQ. 2) X I =AN ( L ) +AN { 2 ) 

IF (NOPT.EQ. 3) X I =A N { 3 ) 


RETURN 

END 


SUBROUTINE GLQU (AUX,ANS) 

THIS SUBROUTINE SETS UP THE L AN CUE RRE-GAUS S INTEGRATION 
IMPLICIT REAL# 8 (A-H,0-Z) 

COMMON /88B/X ( 44 ) , W ( 44 ) 

200 ANS = 0 . 0 DO 
S= 0. ODO 
DO 201 1=1,44 
Y=X ( I ) 

CALL AUX ( Y , Z ) 

20 L S=S + Z#W( I ) 

ANS = S 
RETURN ' 

END 


SUBROUTINE AUX (PHI, FX ) 

THIS SUBROUTINE CALCULATES THE INTEGRAND OF PL, P2, £ P3 
WHERE THE INDUCTION FACTOR = ( ZT-ZT P ) # I NTGRA L OF (PI, P2, CP 3 ) 
NOTE THAT THE INTEGRAND OF P IS GIVEN BY A/B BELOW 
IMPLICIT R E AL #8 (A-H,0-Z) 

COMHON/CCC/ TT, TTP , PI, EL, EP, NR 
COMMON /A AA/ KK, NOPT I 
/ IF tNGPT.EQ.3) PHI=PHl+EP 

(TTrTT r’J PH I K=PH I + 2 . OP 0#P I# ( NB-KK ) / Nd 

A=(tT*# 2}* ( TT-TTP#DCQS ( PHI K) ) + ( TT # ( PH I*DS IN ( PH IiO +DCOS ( PH IK ) )-TTP) 
L *EL##2 

8= ( D SQR T { TT##2 + TTP##2-2.0D0#TT#TTP#DC0S ( PH IK ) +PH I ##2#E L##2 ) ) ##3 
O' 1 #DSQRT( EL##2+<^>#2 ) 

C = A / B LT-f-rp 

FX=DEXP( PHI )#C 
RETURN 
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APPENDIX B 


FORMER APPROACH 

TO THE PROPELLER NOISE PROBLEM 



by C. 

. J Wo an 


LIST OF 

MAIN SYMBOLS 

a 

0 

speed of sound 

t 

observation time 

I 

acoustic intensity 

V 

the velocity of source 

£(r) 

thrust per unit length 

? f 

forward velocity of propeller 

M = 

V/a 0 

X 

field point 

M f = 

V a o 

x = 
0 

: (Xqj + Xq 2 + x^) the distance 




between the field point and the center 

M t = 

r ft/a 
t o 


of the propeller at t 

Mr = 

M-R/R 

Y 

source point 



0 

the instantaneous angle the source 

P' 

acoustic pressure 


makes with y 2 - axis 

P i 

acoustic pressure due to 

8 

rv 

initial angular displacement 


thrust 





f 

tanf = X 02/ X 01 

R = 

X - Y 

T 

retarded time 

X t 

radius of propeller 

ft 

angular velocity of propeller 

T 

total thrust 





[f(t)] = f(T) 



Source Point 



Field Point 


Figure 1 Coordinate System 
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2. Acoustic Pressure Perturbation 


Following Ffowcs Williams and Hawkings, the acoustic pressure perturbation 
at the field point is given by 

P* 


P =" 


l 3 


4TC 3X, ‘ j f R| I - Mr I 


] dS 


th 


where dS is the blade surface area and P is the i component of the 
force acting on air by the blade surface 


3 Noise Due To Thrust 

Let £(r) be the spanwise distribution of the thrust along the lifting line 
per unit length Then 


■ I? 


p =- 


i a 


4tc 3%*. 


~ ;L f g P | R, R —H 
la,R. 2 (l-Mp) a ( st ‘ i- m r 




r — 


R C I- Mg) J 


(*.-&) l-M- 

R I - Mp. 


- P - M ‘]J) 


dr 


where P^ = -£( r) 


P 2 = 0 


P 3 = ° 


Therefore the pressure perturbation due to thrust is 


P t 


where 


t r 1 

4 ^ J ft c r " 9.) 1 f d r 

t 0 


r* = x 

r r t > 


4-0 

^ =-2T 


R = — 


4- * - t 




9C 0 


Tv 
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t CV*) =l*(r*) r* c i - r*) 

1 ft 

T xs the total thrust 
and 


f t ( r > e °) ~ N t r (f-r ) I r*M t xt sint cost — ; s ? 83 

tel 0 -m b ] 3 


4 V Mt Stives 


COS [&] 


[RT(!-M e ] 


3 + r*M t x o Me ~~ n -'f — J ln M — 


- cos f 


CR-’-J (I-MJ 


) 


where the angle 0 is uniformly distributed over (0,2irl Then the proba- 
o 

bility density function for the random variable 8 q is 


■fm = 


i 

2-m 


O ^ 0 < 2JL 


0 


otherwise 


Therefore, the mean square of P 


i 

t 


is 


' 2 

ecu ; 



The intensity 
T = 


Then define 


liKl 

Pq 


4 I 
I - -f2 


FC r, j* r*) Jcr, ) ar drP 


O 0 



where 


2 7t 

FCrfr,*? = j f t (rte») f t or* e„) de„ 


The total thrust is given by 


t = T \ jeer ) r c i- r ) dr 


Letting 


T - i = f iV*) r*(i-r *) dr 

J n 


CD 


the necessary condition that I* be an extremum, subject to the constraint 
J = const (i.e. constant total thrust) is 


T(r*r, ) i c IT ) an =-T r C'~ r > 

J o 

where b is a Lagrangian multiplier 

and Z* (r*) and b can be obtained by solving equations (1) and (2) 
Numerical Method 

The equations needed to solve this problem are 


l - * 


ca> j'iCrVci- r*) dr 


.* 


(B) 


(C) 


CD)' 


= -#rVr*> 

" / 0 ^ 

( 27L 

B (l b r - > =TI^j 'ftO r ffc)4(ir% o)dfio 


i t ( r = Mt r *( i ■ r*) [ r*M t cos * 


CiSj&l 


C2) 




+rM ^ + r ^ Taw ‘ 


-M r ] 
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2 . x 

(E). ££*] = [ (%t cos-^t M f CiA) -t oct S >n 2 -f •+ r^Mt 

^ y z 

-2*xt r^Mt s»n.^ cosC © e + ±*- CRJ) J 
(?) re] = e D + t* - Ift*] 

[GV [l-MpJ = ( I - h| ) - cos^ - r^Mt SirvY SmCSj) 


The solution proceeds as follows 

Given Mp, x*, and T, find b and £*(r*) 


1 The integrals m equations (A) , (B) , and (C) are calculated 
by Gaussian integration (Gauss -Legendre formulas or Gauss- 
Chebyshev formulas) . 


Equation (E) is solved by the Newton Method 

A- I I It 


[^J 


n 


crj„- 


■Kirt,,-,) 






f (Crt) -i TtZS'nV + r Mt 1 

- 2-K-t S<nW / Cos (do -ft * ~ Cfcl) J 

After obtaining [R*] from (2), the value of [R*] is substituted 

into equations (A) , (B) , (C) , (D), (F) , and (G) 


4. After applying Gaussian integration to equation (B) and assign- 
ing n different values for r* (o$r^l, n = the number of points 
used m Gaussian integration), equation (B) can be reduced to 

a system of n equations with unknowns Z*(r*) (i=l,2,--,n) 

til 

[r* are the zeros of the n polynomial used m the Gaussian 
integration] 
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In matrix form 
AL = C 

A. nxn matrix 
C* column matrix 

* ^ T 

L = ii*Cr*), •••, Icr/;) 

5 The Lagrangian multiplier is obtained by substituting L* (r 
obtained ffom (4) into equation (A) . 


* H 
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APPENDIX C 

PRESENT APPROACH TO THE PROPELLER NOISE PROBLEM 
by C. J. Woan 


LIST OF 

A Fourier coefficients of T 

i 

a Q speed of sound 
B number of propeller blades 

C complex magnitude of the 

m th , 

m harmonic 

dD drag per unit length 

dL lift per unit length 

I acoustic intensity 

J Bessel function of order m 

m(2) 

and argument 2 
M = V/a 

o 

Mp = M • R/R 

P-P Q acoustic pressure 

P m root-mean- square pressure 

R = X - Y 

_ 2 2 2 

R = x +y +z 
o J 

r t radius of propeller 


SYMBOLS 

total thrust 
t 

t observation time 

V the velocity of source 

Vp forward velocity of propeller 

W. induced velocity 

X field point 

Y source point 

induced angle of attack 

tj) initial angular displacement 

0 the instantaneous angle the 

source makes with the y^ - axis 

P q fluid density 

T spanwise circulation distribution 

x retarded time 

SI angular velocity of propeller 

[f(t)J = .f(T) 




C-3 


2. Sound Pressure Due to Moving Lifting Line 

The force acting on the air at Y = fix+r, by the moving lifting line, is 

developed as follows 

Consider a blade element at location r 



(dL Cosy - dOs, 71 i') , dF e =< dLsmt + dDcosy 


The force, dF, acting on the air at Y = f2x + r is 

dF - dF^ Z — d p 0 sraQ J + dF e cos 9 'k 


The total sound pressure generated by the moving lifting line is given 
as 


(P~Po)(X,i) 


= _j _3_f r d F<, 

‘ L RC I - Mp.) 


] 


3 Periodicity of the Acoustic Pressure Disturbance 

Consider a propeller whose hub is moving with constant velocity Vp 

Choose fixed axes 0 such that the origin coincides with the hub 

■X.V z 

(i) at time t = 0 

(n) the x axis is parallel to the axis of rotation 
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( 111 ) the hub velocity vector lies in the plane and m x-direction. 

Let the field point "x be at (x,y,z) and a typical source point be situated 
at Y(t) This can be expressed as the vector sum of the hub location and 
the rotational position 


Y<r) = Vp r r (-0 , r(T) = ( 0, rose, r s, ae). 

Here 0 is the instantaneous angle the source makes with y axis This in- 
creases linearly with time, and can be expressed as S2t+<f> where Q is the 
angular velocity of the source, and 4 is the initial angular displacement. 


r = -t - 


CR-J 

CLo 


If attention is confined to source times near t= 0, [R] can be approximated 

by y 

CR] (x-Vp f + ( # - rcos[ejf •+ ( z - rsin[e)) 2 j 


R 0 (i+£) 




r 2 ,.2 , s/%- 

where R 0 — ( X + ^ t 2 ) 


and 


a = £ r* - -2 rC ycoscej 4 - 2 s>n[&] + V f z x -2. % V h zj / Rl 


Therefore, 


tPJ = - -Jr V F r - Y ( + Hs.aCOJ) 


iit - iiy =/2rc l ~M or ) -M r (^cosM + zs,n[«j) 


and 
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where 


M 


or 


X Vf 

O .0 R,o 


and 


M, 


£r 

0~o 


This equation shows that as 0 increases by 2 it, Ut only increases by 2ir 
(1-M ) Hence the basic observed frequency is w=fi/ (l“M or ) 

and therefore P-P o has a period 2 tt/oj 


4 Fourier Analysis 


2ir 


(P-P Q ) (X,t) has a period — , where w = S2/ (r-M Qr ) 


and M 


V x 
F 


or Vo 


th 


Then, defining the complex matnitude of the m harmonic m the usual way 


Cm ( X ) — + <■ b -m ~ 


co 

7t 


IQ 


f M _ Im tot , , 

\ (P-R.X x,t)e dt 

-Jo 


=— L_ ( JL t 
4 TL j axV^ 


where 


I- = — 

TC 


H 


dFt 


PXl-M*) 


J £ dt 


changing variables back to the retarded time, t, gives 


m 


*T- , U) 

” 7C 


j ('JSlj -h ) c]r 


where dt - [ I - M P J dr 

[£] = Ro - a oMor^ -f o a cos ZB] +2 S.ncfl3) 

Since the far field approximation is m force, the only significant 
dependence of F upon x is through the R q term m the exponent. * 
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Therefore 


■ _ jW 1 C r -r n 4*nto(r+-^j dr 

a 0 7c ) L d u * 9^ J 6 ' o 


HtnwC 


{ m to < >n cf> 


C-V it r 7 

m “ ’"T ~ (“'O 

2 k a, fco ( i- M or ) 


( X dF » + 


I- H &£ ft 

M-t r 


-r / M t p \ 

T « (m R;rTU)T) • 

dF 0 . 


Extension to B Blades 

For an assembly of B identical blades , the usual phase arguments show that 
only harmonics which are multiples of B survive, and all the rest cancel 
Hence only harmonics of the propeller blade passing frequency, Bw are 
present m the acoustic field of the complete propeller, and the magnitude 


of the m harmonic is BC : 

mB 


Cm = 


f?o 

■Z •‘■”*“>-*7 

l m 8 co £ 

Zt CLo (l-Mor) 




C"R dF * 


-T d F e ) 


Then, the r m.s pressure is given as 


T v - ?.. K" 


The total intensity at -X is. 


I 
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and the total thrust is 


Tt =5 (r-ft + p o p dr 

0 

Method of Solution 

The problem to solved can be stated as follows 

Find the circulation distribution along the lifting line such that I 

minimum and T constant, 
t 

The solution then proceeds as follows. 

Expand T as Fourier series with coefficients A , then 

3- — X C f \ , , Az , /A 2 • * -) =■ X ( At ) 

In - Tt C A 1 5 A 2 , A 3 • - •) ~ Tt ( At) 


Let J = I - b(T t - 

where b is a Lagrangian multiplier. 

A necessary condition that I be extremum is 


HL = o 


i - I ? -2-j 3 


Solve equation (A) and (B) for A^ and b 



